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Abstract 

We compute the eighth-order fermionic corrections involving two and three 
closed massless fermion loops to the anomalous magnetic moment of the muon. 
The required four-loop on-shell integrals are classified and explicit analytical results 
for the master integrals are presented. As further applications we compute the corre- 
sponding four-loop QCD corrections to the mass and wave function renormalization 
constants for a massive quark in the on-shell scheme. 
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1 Introduction 



In the last about ten years several groups have been active in computing four-loop cor- 
rections to various physical quantities. Among them are the order corrections to the 
R ratio and the Higgs decay into bottom quarks P-[3], four- loop corrections to moments 
of the photon polarization function [IHH] which lead to precise results for the charm and 
bottom quark masses (see, e.g., Ref. [9]), and the free energy density of QCD at high 
temperatures [10]. The integrals involved in such calculations are either four-loop mass- 
less two-point functions or four-loop vacuum integrals with one non- vanishing mass scale. 
In this paper we take the first steps towards the systematic study of a further class of 
four-loop single-scale integrals, the so-called on-shell integrals where in the loop massless 
and massive propagators may be present and the only external momentum is on the mass 
shell. 

On-shell integrals enter a variety of physical quantities, where the anomalous magnetic 
moments and on-shell counterterms are prominent examples. The first systematic study of 
two- loop on-shell integrals needed for the evaluation of the on-shell mass and wave function 
renormalization constants {Z^^ and ^°^) for a heavy quark in QCD has been performed 
in Refs. [II1I12]- Already a few years later, in 1996 the analytical three-loop corrections 
to the lepton anomalous magnetic moment ai became available [13] . This result has been 
checked in Refs. [T^UTS] . In Refs. [T^|[T6] the three- loop on-shell integrals have been applied 
to QCD, namely the evaluation of and 2^^^. The calculation of Ref. [H] has confirmed 
the numerical result of [17l[T8] which has been available before. Both and have 
also been computed in Ref. [I5]. Further application of three- loop on-shell integrals are 
discussed in Refs. [T91I20]. There is no systematic study of four-loop on-shell integrals 
available in the literature. Nevertheless, some four-loop results to the anomalous magnetic 
moment of the muon, a^, have been computed analytically, in particular contributions 
from closed electron loops. E.g., the contribution where the photon propagator of the 
one-loop diagram (see Fig. [1]) is dressed by higher order corrections has been considered 
in several papers [2T1 - I27] . Four-loop corrections where one of the two photon propagators 
of the two- loop diagram is dressed by higher orders has been considered in Ref. [281129] . 
Contributions where both photon propagators get one-loop electron insertions are still 
missing. This gap will be closed in the present work. Let us mention that all four- and 
even five-loop results for a; are available in the literature in numerical form |27 [[5DH55] 
(see also the review articles [MIES])- 

In this paper we take the first step towards the analytical calculation of four-loop on-shell 
integrals by considering the subclass with two or three closed massless fermion loops, 
which are marked by a factor rii. Thus we are concerned with four- loop terms proportional 
to nf and nf which we consider for three physical quantities: the anomalous magnetic 
moment of the muon, a^, the on-shell mass renormalization constant, Z^^, and the on- 
shell wave function renormalization constant, 2'^^, for a massive quark. For the latter 
QCD corrections to the quark two-point functions are computed whereas for the former 
muon-photon vertex diagrams have to be considered. Some sample Feynman diagrams are 
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Figure 1: Sample Feyman diagrams for the photon- muon vertex contributing to a^. Wavy 
and straight hues represent photons and fermions, respectively. In this paper we consider 
the contribution where at least two of the closed loops correspond to massless fermions. 
The last diagram in the second line is a representative of the so-called "light-by-light" 
contribution. 

given in Figs. [1] and O The precise definition of these quantities is provided in Sections [3] 
andU 

The outline of the paper is as follows: in the next section we provide details of the four-loop 
on-shell integrals needed for our calculation. In particular, we identify all master integrals 
and provide analytical results in Appendix A. The renormalization constants and 
are discussed in Section |3] and Section |4] is devoted to the anomalous magnetic moment of 
the muon. We discuss the relation between the MS and on-shell fine structure constant 
and provide analytical results for a^. Finally, we conclude in Section |5l Appendix B 
contains the analytic results for the relation between the fine structure constant defined 
in the MS and on-shell scheme. 
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Figure 2: Sample Feynman diagrams for the QCD corrections to the fermion propagator 
contributing to and 2'^^. Curly and straight lines represent gluons and fermions, 
respectively. In this paper we consider the contribution where at least two of the closed 
loops correspond to massless fermions. 

2 Four-loop on-shell integrals 

In this Section we present the setup used for the calculation and discuss the families of 
four-loop on-shell integrals needed for the nf and nf corrections for 2'°^, and a^. 
Since all three cases reduce to the calculation of corrections to the fermion propagator we 
consider in this Section the corresponding two-point function. 

After the generation of the diagrams with QGRAF we use q2e [37l|38] to translate the 
output into a FORM [39] readable form. In a next step exp [371 138] is applied to map 
the momenta to one of five families. During the evaluation of the FORM code we apply 
projectors and take traces to end up with integrals which only contain scalar products in 
the numerator and quadratic denominators. 

In the next step we have to reduce all occurring integrals to a minimal set of master 
integrals. This is done using two different programs in order to have a cross check for 
the calculation. On the one hand we use crusher [IQ] and on the other hand the C++ 
version of FIReE Both programs implements Laporta's algorithm |42] for the solution of 
integration-by-parts identities [13]. We find complete agreement for the expressions where 
the physical quantities are expressed in terms of master integrals. 

Let us mention that we have performed our calculations for general gauge parameter 
which drops out once the four- loop results for and are expressed in terms of 

^The Mathematica version of FIRE is publicly available [4T| . 



4 



master integrals^] 

Altogether we end up with 13 master integrals. Seven of them (shown in Fig. [3]) are 
products of one- and two-loop integrals whereas the remaining six integrals (cf. Fig. |4]) 
request a dedicated investigation. We calculate them using the Dimensional Recurrence 
and Analyticity (DRA) method introduced in [44]. In order to fix the position and order 
of the poles of the integrals, we use FIESTA [l5lll6]. The remaining constants are fixed 
using the Mellin-Barnes technique [TFl - IST] . In order to express the results in terms of 
the conventional multiple zeta values we apply the PSLQ algorithm [52] on high-precision 
numerical results (with several hundreds of decimal digits) o 

The analytic results for the integrals in Fig. [Hare listed in Appendix A. Results in terms 
of Gamma functions for the integrals in Fig. |3] are easily obtained recursively using the 
formulae from the Appendix of Ref. [19]. For convenience also these results are given in 
Appendix A. 

All results have been cross-checked numerically with the help of FIESTA [46] where an 
accuracy of at least four digits has been achieved. 



3 Fermionic nf and contributions to Z^^ and 

Both and are obtained from the fermion two-point functions S(g) which can be 
cast in the form 

E(g,mq) = mgSi(g^,mg) + (j^ -mg) E2(g^,mg) . (1) 

Here rriq represents a generic quark mass whereas bare, on-shell and MS quark masses are 
denoted by m°, Mg and rhg. 

The derivation of ready-to-use formulae for Z^^ and is discussed at length in Refs. [HI 
[T5] . Thus, let us for convenience only repeat the final formulae which are applied in our 
calculations. They read 

Z^^ = 1 + Si(M2,M,), (2) 

{Z?T' = l + 2M;^Sl(g^M,)_^^^ + S2(M^^M,). (3) 



The expressions on the right-hand side are computed by introducing the momentum Q 



^Note that and have to be independent of the QCD gauge parameter f whereas we expect 
that the nj and n/-independent terms of do depend on ^. 

■^Let us mention that the numerical evaluation of the factorizable four-loop master integrals for a; 
which reduce to the evaluation of the corresponding three-loop master integrals in higher orders of e was 
undertaken in Ref. |53j as a warm-up before a future full four-loop calculation. This was done with the 
method of [42] based on difference equations. The achieved accuracy of several dozen of decimal digits 
was not enough for using PSLQ. 
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Figure 3: Master integrals for the and nf contribution which are easily obtained by 
applying one- and two- loop formulae, see e.g., Ref. [H]. Solid lines carry the mass M 
and dashed lines are massless. For Li to Lq we have = where q is the external 
momentum; Lj is a vacuum integral. 




Figure 4: Non-trivial master integrals contributing to the nf contribution. The same 
notation as in Fig. [3] has been used. 
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with = via q = Q{1 + 1) which leads to the equation 



4M2 

. (4) 

Hence, to obtain Z^^ one only needs to calculate Si for = M|. To calculate Z^^, one 
has to compute the first derivative of the self-energy diagrams. Note that the renormal- 
ization of the quark mass is taken into account iteratively by explicitly calculating the 
corresponding counterterm diagrams. 

We write the perturbative expansion for in terms of the renormalized strong coupling 
as (7^; is the Euler-Mascheroni number) 



IT \ Ait J " \ 7T J \A7T 



This allows us to take the ratio between the on-shell and MS [5H - I56] mass renormalization 
constant which is given by 



mg(/i) 



^os 



+ O {al) (6) 
The coefficients 6zm are by construction finite. 

In the case of we choose the bare coupling as expansion parameter which in many 
applications turns out to be convenient. Furthermore, the dependence on fi/Mg can be 
written in factorized form which leads to shorter expressions. Thus we have 

where each term ^Zg"^ contains a factor (yU^/Mg)"*^. 

We refrain from repeating the one-, two- and three-loop results for Z^^ and since 
analytical expressions for general colour coefficients are available in the literature [TIHI6]. 
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We split the four-loop coefficient according to the number of closed massless fermion loops 
and write {i G {m, 2}) 



with an analog notation for 6zm . 

In th 

read 



In the following we present analytical results for 6zm'^\ 6zm'^\ bz'^'^'^ and bZ^^"^^ which 
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where Im = In Cn is Riemann's zeta function, ai = In 2 and a„ = Li„(l/2) [n > 1). 

In the case of QCD the colour factors take the values Ca = 3, Cp = 4/3 and T = 1/2. In 
Eqs. f lTOj) and f lT2|) the contributions from closed heavy quark loops are marked by Uh = 1 
which has been introduced for illustration. 

In order to get an impression of the numerical size of the newly calculated terms we 
evaluate for /i = Mg. After inserting the numerical values for the colour factors we 
obtain {As = as(Mq)/vr) 

z^^ = 1-A,1.333 + A^ (-14.229 - 0.104n/, + 1.041 n;) 

+ Al (-197.816 - 0.827 Uh - 0.064 + 26.946 - 0.022 n^^n; - 0.653 n^) 
+ Af {-A3 A65 -0.017 UhU^ + 0.678 + ...) +0 (Al) , (13) 

where the ellipses indicate rii independent contributions and terms proportional to ni 
which have not been computed. One observes that the nf contribution at two loops and 
the nf contribution at three loops are quite small. This is in contrast to the linear ni 
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terms which can become quite sizeable. E.g., setting rii = 5, which corresponds to the 
case of the top quark, we obtain (for 71^ = 1) 



z^^ = 1-^,1.333 + ^2 (-14.332 + 5.207„J 

-198.707 + 134.619„, - 16.317„2 



+ Al (- 

+ At (-1087.060^2 + 84.768„3 + ...)+ O (Al) . (14) 

At two-loop order the rii contribution is only a factor of three smaller than the ni- 
independent term, however, with an opposite sign. At three loops the linear-ri; term 
has almost the same order of magnitude than the constant contribution but again a dif- 
ferent sign. It is remarkable that for rii = 5 the coefficient of the four-loop nf term is 
more than a factor of five larger than the nrindependent term at order a^. 



Let us finally compare our results with the approximate expressions obtained in Ref. 
in the large-/3o approximation. In Ref. [57] one finds for the quantity Mg/rhq{rhg) the 
result (a, = as{rhg)/7i) 



mq{mq) 



large - 



1 + a,1.333 + al (17.186 - 1.041^;) 
+ al (177.695 - 21.539n, + 0.653^2) 

+ at (3046.294 - 553.872^; + 33.568nf - 0.678nf) , (15) 



where for the renormalization scale fi = ruq has been chosen. The coefficients of Eq. (fTS!) 
should be compared with our findings which read 

= 1 + a,1.333 + a^ (13.443 - 1.041n/) 



mq[mq) 



+ al (190.595 - 26.655n; + 0.653nf) 
+ at (co + ciui + 43.396^2 - 0.678nf) , (16) 

where Cq and Ci are not yet known. By construction one finds agreement for the coefficient 
of nf since it has been used as input in Ref. ^7\. As far as the nf term is concerned the 
exact coefficient is predicted with an accuracy of about 30%. 



4 Fermionic and nf contributions to 

It is convenient to introduce the form factors Fi and F2 of the photon-lepton vertex as 

T''{q,p) = F,{q')r + ^^^a,.q\ (17) 
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where q is the incoming momentum in the photon hne and Mi is the lepton mass. The 
anomalous magnetic moment is given by 



In Eq. f lT7|) also the momentum p = [pi +p2)/2 has been introduced where Pi = P2 = 
are the momenta flowing through the external fermion lines (see Fig. [1] for the directions 
of the momenta). 

The evaluation of a/ requires that r^(g,p) is computed in the limit g — ?■ 0. Due to the 
factor q'^ in front of F2 in Eq. f fT7|) one has to perform an expansion of r^{q,p) up to linear 
terms in q which can be written as 

r^(g,p) = X'^ip) + q,Y'^-'ip) + 0{q') , (19) 

with = M^. F2 is conveniently obtained after the application of a projector given by 
(see, e.g., Ref. [58] ) 



1 ^ 

2Mf{D -1){D -2) ^ 



2 

(20) 



and thus ai is reduced to the evaluation of on-shell two-point functions as described in 
Section [2l 

We define the loop expansion of ai in analogy to Eq. ([5]) (with as replaced by the fine 
structure constant) and introduce the same splitting according to the number of massless 
lepton loops as in Eq. (jH]). 

The Feynman diagrams contributing to ai respectively the coefficients of a" and raf can 
be subdivided to two classes: (i) the one where the external photon couples to the lepton 
at hand and (ii) the one where it couples to a lepton present in a closed loop. Sample 
diagrams are given in Fig. [1] In the following we refer to the diagrams of class (ii) as 
"light-by-light" contribution in analogy to the corresponding hadronic part. 

In this paper four-loop corrections contributing to class (i) are evaluated which contain 
two or three closed massless fermion loops. They are used in order to compute electron 
loop contributions to neglecting terms of order M^/M^. 

For the diagrams in class (i) we can proceed as follows: In a first step we renormalize 
the fine structure constant in the MS scheme, a(/i). The corresponding renormalization 
constant is easily obtained from the one for after specifying the colour factors to QED. 
The MS scheme has the advantage that the electron mass can be set to zero (which is 
not the case for the diagrams in class (ii)). After renormalizing the muon mass in the 
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on-shell scheme we obtain a finite expression for whicli sliows an explicit dependence 
on In (/iVM^). 

In a next step we replace by its on-shell counterpart using the corresponding relation 
up to three loops. It can best be calculated by considering the photon two point function 



I / dxe"'%Q\f{x)f{x)\Q) 



(21) 



and employing the on-shell renormalization condition n(g^ = 0) = 0. The form of the 
renormalization condition reduces the problem to the calculation of two-scale vacuum 
integrals at three loops. Note, that for the renormalization of the fermion masses in the 
on-shell scheme the dependence on both masses has to be taken into account. In the limit 



■C we obtain 



(see also Refs. [25l|27l|59]) 
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with Lfj_ = In(/i^/M^) and Le = In(/i^/Mg). The ellipses in the coefficient of (a/n)'^ 
indicate terms which we left out since they are irrelevant for the contribution discussed 
in this paper. The complete result containing the exact dependence on M^/M^j, is presented 
in Appendix B. Note that the result in Eq. ( 122|) can be obtained from the one provided 
in Ref. [27] where the relation between a(/i) and a is given for one massive lepton. 

Also in the case of ai we refrain from listing the lower-order results which can be found 
in the literature [T3 | 1321435] . Rather we concentrate on the new correction terms at four 
loops. Adopting the notation from Eq. (jH]) we obtain the following results for the nf 
contribution 
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(23) 



where L^e = ^^i^'^/^e)- The approximate results have been obtained with the help 
of [60] M^/Me = 206.7682843(52). The result in Eq. ([23]) agrees with the one in Ref. [281 



In the case of the nf contribution we split a^f^"* into two parts. The first one 
corresponds to the diagrams containing two closed fermion loops and the second one 
originates from diagrams with three closed fermion loops where one of them is a 
muon and two are electron loops. Thus, we have 
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with 
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(25) 



,(42)6 



agrees with Ref. [28|[29]. Analytical results for al^^^^" are not present in the literature 
since corrections originating from diagrams as the third one in the first row of Fig. [1] have 
not been considered yet. However, we can perform a numerical comparison with the 
results from Refs. [30ll33li which reads 



,(42)a 



-3.642 04(1 12) 



(26) 



There is a good agreement with the analytic result in Eq. ( 1241) . The deviation can be 
explained by corrections of order MJM^ ^ 0.005 or (Me/M^)^ In^ M^/M^ 
which are absent in our analytic expressions. 



5 Conclusions 

In this paper the first steps towards the evaluation of four-loop on-shell integrals have 
been undertaken. As an application within QCD we have computed the contributions 
involving two massless quark loops to the on-shell renormalization constants and 
Z^. As an application in QED we have considered the contribution from four- loop 
diagrams involving two or three closed electron loops to the anomalous magnetic moment 
of the muon excluding, however, the light-by-light contribution. 

We describe in some detail the techniques and the programs which have been used for the 
calculation. We are confident that they are generic enough to be applied to the n\ and 
non-fermionic contribution. The only bottleneck might be the analytic evaluation of the 
master integrals so that maybe numerical methods have to be applied. 

^In in Ref. |33| the contributions from closed electron and muon loops are always added whereas in 
our result at least two closed electron loops are present. We are deeply grateful to the authors of Ref. [33] 
for providing us the results for the contributions containing only electron loops Eq. p6|) . 
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Appendix A: Analytic results for the master integrals 

In this appendix we provide the analytic results of all master integrals where we assume 
an integration measure d^k / {iTi)^/"^ with D = A — 2e. Furthermore we write scalar 
propagators of particles with mass M in the form l/{-k'^ + M^). For convenience we set 
M = 1 in the final result since the dependence on M can easily be restored. 

The analytic results for the integrals in Fig. [3] read 
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r(D - 2)2r(2D - 5) 

,2^ /n ,x4. 



^ _ r(3-D)2r(2-f)'r(f -i)*r(2D-5)2 
r(D-2)2r(if -3)' 

_ r(5 - 2D)T (4 - ^) r (f - 1)^ r(4D - 9) 

" r(2D-4)r(^-5) ' 

_ r(6 - 2D)v (5 - f ) r (2 - f )' r (f - 1)' r - 4) t{ad - 11) 

^ ~ r(4 - D)T{D - 2fT{2D - 5)1 - 6) 

r(6 - 2L))r(3 - D)T (2 - f ) r (f - 1)^ T{AD - 11) 

r(D-2)r(f -3)r(f -6) ' 

r(7 - 2D)T (2 - f )^ r (f - lfT{AD - 13) 

r(Z} - 2)3r (f - 7) ' 

^ r(6 - 2D)v (5 - f )^ r (2 - f )^ r (f - 1)^ r (f - 4) 
^ r(io-3D)r(i:)-2)2r(f) 

The analytic results for the integrals in Fig. |4] read 
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(32) 



5 25 /205 177r2\ ^ f 323 8571^ 79(3 \ _i [ 55241 



24e4 24e3 \QQ 72 J \ 72 18 J \ 1152 

^ 4097r2 ^ 395C3 ^ vr^^^o 733351 ^ 41997r2 1943(3 57r^ 4997r2C3 



288 18 8 / V 3456 288 72 8 54 
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_ 407C5\ 1 / 14346449 3830457r2 19057(3 4377r4 249571^(3 2035(5 
Q~)^ " V 41472 ^ 3456 ^ 72 ^ ~96 54 6 

20 2 77r6 2285C|\ ^ ( 391938053 35179637r2 1751323(3 9 3 3 477r^ 

^ 11340 '2n~)^ ~\ 124416 ^ 10368 ^ 864 ^ 1440 

_ 4437r2C3 21905C5 20277r6 11425C| 24777r^C3 9223712^5 11681C7A ; 
216 ^ 24 ^ 2268 27 90 90 ^ 42 y 

+ O {e') , (33) 



23 /433 297r2\ ( 297 19l7r2 275(3 \ _i ( 22765 



e^^^^Ms = — - ^ - — + ^ 



12e4 8e3 V 48 36 7 V 32 24 18 y V 64 

, 71777r2 2 ^ 2273<3 ^ 1257rn / £1105 SOSSvr^ 
H 247r ai H ^ ^ — 1 3247r ai 



144 12 72 y V 128 32 

, 105463(3 , 27477r^ , on 2 2 , ,n 4 , 15957r^C3 , 3223(5 , \i 
+ + + "1 + + ^4 + + )^ 

16944559 , 21673l7r2 ^ , 146091(3 , 437577r^ , ^^^^ ^ ^ 

27067r oi H -— 1 — — h 10807r 



768 192 16 1440 

2 3 on 5 , 9785712(3 52351(5 , x.^o... , ...^.S3 



167r% 800 2 3 on 5 ^ 97857r2(3 52351(5 , 1123397r6 15125( 

+ 540ai H TT a? — 80ai H 

^ 3 3 ^ ^ 36 20 22680 54 

, , , 6869772 1 5898057r2 „ 4851365(3 

+ 12960a4 + 9600a5 — — + — — 1816571^01 ' ^'^ 



).--(; 



512 128 96 



9485371^ . . 44 93 5 3364157rV3 

+ 902071^0? + 4510at + 727r^ai - 36007r2a? - 1080a? + ^ 

960 111 11 216 

8849321(5 , 3515997r6 . „ „ , 400a? , . > 
+ 32640^6 - - 1047r4a? + 7447r2at + + 9447r2ai(3 

375097(2 687571^(3 934677r2(5 652775(7 

^ + ^ + — + — - + 108240a4 + 129600a5 + 96000a6 

36 108 90 42 

+ 18567rV)e^ + (e^) , (34) 



. _ 1 7 /lO 137r2\ 2 / 61 737r2 118(3\ 1 / 851 837r2 



6e4 66-^ V 3 18 y V 6 18 9 y V 4 1 

^ 637(3 ^ 377r^\o f 1486 1 34677r2 ^ 1003(3 ^ 112l7r^ , 189471^(3 



+ 



10 y V 8 36 18 60 27 

16018(5 \ 1 / 613975 2598171^ 68293(3 837r4 95597r2(3 79891(5 



15 y V 48 24 36 24 27 15 

^ 5950l7r6 ^ 17704(|\ 2 _ 7539347 3823497r2 482627(3 4266597r^ 



2835 27 y V 96 48 24 720 

37577r2(3 2525(5 4320l7r6 88585(| 1720471-^(3 20643471^(5 
^ 54 ^ 6 ^ 420 ^ 27 ^ 45 ^ 45 



15 



+ H5|!?Ci),3 + o(.^). (35) 



,*.«M, = -±-l^-('l^*^\-^-(3+l^ + ^e-^ - (-2,0 ■ '82='' 



3e4 2e3 V 6 9 y V 3 9 J \ 36 

^.2 ^ 464C3 ^ 197rn f 5091 , 281l7r=^ , , 14797(3 , ITn^ 

- 327r ai H e — 4327r ai H 

3 45 ; V2 8 9 90 

320 2 2 leOaf 33271^(3 1556(5 \ i f 55049 956937r2 

- 3msA. + HHHCi _ }^ + i440.'a; + 720a? + H2l!^ - 5H22,=„3 

2 45 9 9 

320a? 10467r2C3 , 7727r6 2648(1 \ o 

- + ^ - 8246(5 + + + 17280a4 + 12800a5 J 

45814 1 3294677r2 „ 938425(3 99797r^ 36080 o 2 
- 242207r ai + — -— + -^—ir 



4 32 12 36 3 

+ + 240.'a, - 4800.=a? - 1440a5 + _ + 44800.. 

769047r6 1 7« 4 2 ^ a^. 2 4 ^ leOOof ^ 4160 2 155668(| , 43047r^(3 

"^4^"^^^'' ai + 9767r a, + + ^vr ai(3 ^ + 135 

+ !^ + 1^ + 144320., + I72800a5 + 128000ae + 

+ O (6^) , (36) 

^4.^,^ ^ 1 11 _ /l^ , _ + 1^ , / 28699 



12e4 24e3 V 16 36 7 V 96 72 9 7 V 192 

^ 657r2 ^ 559(3 ^ 1497r^\ _ / 144429 _ 147557r2 ^ 227(3 , 1937?!^ 



16 9 90 y V 128 288 2 180 

^ 11187r2(3 ^ 7604(5 \i ( 5327075 3 730877r2 45889(3 , 7497r^ 



27 15 y V 768 576 36 40 

^ 72677r2(3 ^ 49426(5 ^ 140537r6 ^ 14063(|\ 2 58275695 6258597r2 



27 15 1620 27 y V 1536 128 

1192693(3 1681437r^ 295171^(3 1826897r6 182819(| 

h h 58O5C5 H V 

72 720 6 ^ 3240 54 

510137r4(3 , 988527r2(5 ^ 102171 1(7\ 3 



270 45 42 ' ^ ^ ' ^ ' 



1 17 /149 137r2\ _2 / 433 1497r2 23(3 \ _i / 3817 



e - — + V — e - 



12e4 24e3 V 48 36 / V 96 72 9 / V 64 
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^ 52l7r2 ^ 173C3 ^ 34l7r^\ f 97165 941971^ ^ 1367(3 ^ 166771^ 



144 9 180 y V 128 288 18 180 

^ 659712(3^ 5939(5^1 f 4640963 8046l7r2 ^ 1717(3 , 8337r4 ^ 35277r2(3 



+ 



27 15 y V 768 192 36 360 27 

30599(5 5779l7r6 734(|\ ^ ( 61900879 4102837r2 55357(3 



15 5670 27 y V 1536 128 24 

193 9 5 77r^ 84897r2(3 45833(5 110691l7r6 14197(| 2 1 6377r^(3 

720 ^ 54 ^ 30 ^ 22680 ^ 54 ^ 135 
754 077r2(5 52201 7(7 \ 3 



+ ^^^^ + ™IJe3 + 0(e^), (3^ 
with S6 = YZ=i Er=i(-1)"^VK^) = 0.98744 



Appendix B: Relation between and a 



In this Appendix we present the result for the relation between the fine structure constant 
defined in the MS and on-shell renormalization scheme involving two massive leptons with 
masses mi and m2. We label contributions from leptons with mass mi and m2 by Uh and 
ni, respectively. Our result reads 

2 



a{ij) 



a 



1, a 1, , 1,22 h 15 

-hni- + < -hhnhni + -h Ui + -r + T7 
3 7r 9 9 V 4 16 



ni 



3^ 3 



+ ni 



384 







'5/2^ 

. 24 


k (3 

2 S3 


57r2 


77 \ 


32 192 


24 


^ 5767 


79/2 (3a;2 - 




5/1/2 


1152 




24 



79/2 

144 



7(3 7r2 695 
64 ^ y ^ 648 



384 



(-128x^ - 15x2 - 71) ln2( 



X) 



+ \ {-x^ + + X-1) Li2(l -x) + ^{x^ + x^ + X + 1) Li2(-x) 



+ 



{5x^ + 3x^ + 3x2 ^ 



256x3 



1 



+ - (x^ + x^ + X + 1) ln(x) ln(x + 1) + 
3 



(^(Li2(l - x) + Li2(-x)) ln(x) - 2Li3(l - x) - Li3(-x; 

(5x^ + 3x^ + 3x2 ^ Y^2(^^^ ^ 



512x3 



+ 



405x^(3 + 11527r2 (x^ 



X 



5994x2 + 243x(3 - 5126 



10368 



+ <nh ^ ni,mi ^ m2, x f-)- 



with X = mi/m2, /fc = ln(/i2/m^) and Oi = In 2. 



X J 



(39) 
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